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Abstract. The system of subdwarfs G89-14 is one of the most metal-poor multiple stars with an atmospheric metal 
abundance [m/H] = —1.9. Speckle interferometry at the 6-m BTA telescope has revealed that G89-14 consists of 
four components. Measurements of the magnitude difference between the components and published data have 
allowed their masses to be estimated: M A ~ 0.67 M Q , M b w 0.24 M s , M c « 0.33 M©, and M D « 0.22 M Q . The 
ratio of the orbital periods of the subsystems has been obtained, 0.52 yr : 3 000 yr : 650 000 yr (1:5769:1250000), 
indicative of a high degree of hierarchy of G89-14 and its internal dynamical stability. The calculated Galactic 
orbital elements and the low metallicity of the quadruple system suggest that it belongs to the Galactic halo. 



1. INTRODUCTION 

Multiple star systems (with ^ 3 components) are cur- 
rently believed to be a natural result of star formation 
(Kroupa et al. 2003; Delgado-Donate et al. 2003). At least 
8% of the solar-type stars consist of three or more compo- 
nents (Tokovinin 2008). Since multiple systems are char- 
acterized by additional parameters compared to single and 
double stars (the angles between the orbital planes of var- 
ious subsystems, eccentricity ratios, etc.), their study al- 
lows more stringent constraints to be imposed on the star 
formation process. Star formation models must predict the 
properties of such objects, while the theory of dynamical 
evolution must predict their survivability under both in- 
ternal (tidal interactions) and external (the influence of 
the Galaxy's nonuniform gravitational field) effects. 

Investigating metal-poor multiple stars ([Fe/H] < —1) 
is important in understanding the stability of star systems, 
because the stars with a low atmospheric metal abundance 
were generally formed at the epoch when our Galaxy was 
born, i.e., more than 10 Gyr ago. Studying multiple sys- 
tems that have survived over such a long period allows one 
not only to test various criteria for dynamical stability of 
multiple stars but also to impose constraints on the mass 
distribution in the Galaxy. 

In 2006 and 2007, we conducted a speckle interfer- 
ometric survey of 223 nearby F, G, and early-K subd- 
warfs (< 250 pc) with low metallicities ([Fe/H] < —1) 
and large proper motions (p, > 0.26"/yr ) at the 6-m 
BTA telescope of the Special Astrophysical Observatory, 
Russian Academy of Sciences (Rastegaev et al. 2007, 2008; 
Rastegaev 2009). The goal of the survey of Population II 
stars, which was carried out with the diffraction-limited 
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angular resolution of the 6-m telescope (0.023" at 550 
nm), was to expand the database of double and multi- 
ple old stars and to determine the orbital parameters and 
properties of their components. One of the results of this 
survey was the discovery of the quadruple system of sub- 
dwarfs G89-14 (HIP 35756). 

In this paper, we present the fundamental parameters 
of this system obtained from our observations and pub- 
lished data. 



2. OBSERVATIONS 

The speckle observations of G89-14 were performed at 
the 6-m BTA telescope of the Special Astrophysical 
Observatory, Russian Academy of Sciences, in December 
2006 and March 2007 (Rastegaev et al. 2007, 2008). 
In the observations, we used a system (Maksimov et 
al. 2009) based on a 512 x 512-pixel EMCCD (Electron 
Multiplying CCD) with a high quantum effciency and lin- 
earity, which allowed objects with a magnitude difference 
between the components Am < 5 m to be discovered with 
the diffraction-limited resolution of the 6-m telescope. The 
size of the detector field, 4" made it possible to detect the 
secondary components at angular distances as large as 3" 
from the primary star. 

In December, under good weather conditions (seeing 
< 1"), we accumulated 500 20-ms exposures for G89-14 
through a 800/100 nm filter (the first and second num- 
bers give the central wavelength of the filter passband 
and the passband FWHM, respectively). In March, the 
weather conditions were not optimal for speckle observa- 
tions (w 3"). We obtained 2000 speckle images with 20- 
ms exposures through each of the two filters, 550/20 and 
800/100 nm. 



2 Rastegaev: The Most Metal-Poor Quadruple System 

Table 1. Results of speckle measurements for G89-14 
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A/AA, nm 


Epoch 


0.979 


0.009 


0.8 


0.4 


4.1 


0.4 


800/100 


2006.94455 


0.982 


0.005 


0.8 


0.4 


4.3 


0.1 


800/100 


2007.24040 


Unresolved 


550/20 


2007.24040 



We calibrated the measurements based on the so-called 
"standard" pairs — binary systems with well known sepa- 
rations between the components and position angles. The 
technique for determining the relative positions and mag- 
nitude differences of the components of the objects under 
study from speckle interferograms averaged over a series 
of power spectra was described by Balega et al. (2002). 
The accuracy achieved with this technique is 0.02" 1 for 
the magnitude difference, 1 mas for the angular separa- 
tion, and 0.1° for the position angle. 

The results of our measurements for the speckle inter- 
ferometric subsystem of G89-14 are presented in Table 1. 
The March observations through the 550/20 filter did not 
reveal the speckle intcrferometric component at a distance 
of 0.98" from the known spectroscopic pair (see below), 
because its magnitude in this spectral range was fainter 
than that of the SB1 system by more than 5 m . At the 
same time, the observations through the 800/100 filter in 
March 2007 confirmed, within the error limits, the results 
of our intcrferometric observations in December 2006. 



3. G89-14: A SYSTEM OF FOUR SUBDWARFS 

Initially, G89-14 (TYC 763-801-1; HIP 35756; WDS 
07224+0854; NLTT 17770) with coordinates 07 h 22 m 31!5 
+08°49'13.0" (J2000.0) was known as an SBl-type spec- 
troscopic binary system with a period of about 190 days 
(Latham ct al. 1988, 2002). Subsequently, Allen et al. 
(2000) found a companion to this pair with a common 
proper motion located at an angular distance of 34". 
Unfortunately, these authors did not give the position an- 
gle of this companion. In the Aladin database, there is 
only one star at a distance of 34" that we took as the com- 
ponent with a common proper motion (Fig. 1). Allen et 
al. (2000) refuted the physical connection of G89-14 with 
G89-13 (HIP 35750), because, having a common proper 
motion, the parallax and radial velocity of G89-13 dif- 
fered significantly from those of G89-14. In 2006, using 
the BTA telescope, we discovered a speckle interferometric 
component at an angular distance of « 1" from the spec- 
troscopic binary (Rastegaev et al. 2007, 2008). According 
to the Hippacros data, the heliocentric distance of G89- 
14 is w 170 pc (ESA 1997). The distance determined by 
Carney et al. (1994) from the photometric parallax is 94 
pc and is most likely an underestimate, because it dis- 
regards the presence of additional components. At a dis- 
tance of 170 pc from the Sun, the probability that the 




Fig. 1. Image of G89-14 from the DSS2 red survey: the 
triple subsystem consisting of the spectroscopic pair and 
the speckle interferometric component is at the center; 
the component with a common proper motion is located 
34" to the northwest. In the figure taken from the Aladin 
database, the north is at the top and the east is to the 
left. 



speckle interferometric companion at an angular distance 
of 0.98" and with a magnitude difference of about 4™ is 
an optical projection is insignificant. Hence, we concluded 
that the discovered companion and the system G89-14 
are physically connected. Their gravitational connection 
is also confirmed by the repeated interferometric obser- 
vations in March 2007. If the component were an optical 
projection, then, given the large proper motion of G89- 
14, n w 0.3"/yr (Carney et al. 1994; ESA 1997), the rel- 
ative position of the component would change by w 0.1" 
in the time between our observations (from December to 
March) , which exceeds the errors of our speckle measure- 
ments. Thus, four gravitationally bound components of 
the system are known at present (Fig. 2). The metallicity 
of G89-14 is [m/H] = -1.9 (Carney et al. 1994). Analysis 
of the literature has shown that none of the quadruple 
stars known to date has such a low metal abundance. 
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Fig. 2. Schematic representation of the hierarchy levels 
for G89-14: A-B is the spectroscopic subsystem, AB-C is 
the speckle interferometric subsystem, and ABC-D is the 
subsystem with a common proper motion. 



4. ESTIMATES OF THE COMPONENT 
MASSES 

Accurate stellar masses can be determined only in bi- 
nary and multiple systems by constructing their orbits. 
In the remaining cases, only approximate estimates based 
on various theoretical and empirical relations (e.g., the 
mass-luminosity relation) are possible. All orbital param- 
eters are known for none of the subsystems of G89-14. A 
spectroscopic orbit was obtained for the shortest-period 
subsystem with a period of 190.49 ± 0.29 days (Latham 
et al. 2002). However, in this case, the inclination i of the 
orbital plane to the plane of the sky remains unknown. 

Using evolutionary tracks, spectroscopic (Carney et al. 
1994) and photometric (Allen et al. 2000; Monet et al. 
2003) data, and our speckle observations, we estimated 
the masses of the components in the quadruple system 
under consideration. 

Based on the tracks from Vandenberg and Bell (1985) 
and the spectroscopically measured metallicity and tem- 
perature of the primary component (5750 K) of the spec- 
troscopic subsystem of G89-14, Carney et al. (1994) esti- 
mated its mass to be 0.67 M . The mass function that 
they obtained from spectroscopic observations, 



Ml sin 3 



M B y = °- 0173 < 



(1) 



allows the mass of the secondary spectroscopic component 
to be determined to within a factor of sin 3 i. Substituting 
M A = 0.67 M Q into Eq. (1) and setting sin 3 i = 1, we 
obtain the lower mass limit Mg > 0.24 M . 

The mass of the speckle interferometric companion 
was estimated from the total magnitude of the triple 
subsystem (the SB1 pair and the interferometric com- 
ponent), 10.0™ in the / band (Monet et al. 2003), and 
from the magnitude difference between the spectroscopic 
binary and the companion we discovered (Table 1). For 



our calculations, we took Am = 4.2 in the 800/100 filter, 
which roughly corresponds to the / band. The apparent 
magnitude of the speckle interferometric component was 
14.22™. Based on the absolute magnitude of the compan- 
ion (Mi = 8.07 m ) and using the evolutionary tracks for 
[Fe/H] = -2.0 from Baraffe et al. (1997), we estimated the 
mass of the interferometric component to be ~ 0.33 M . 

Allen et al. (2000) gave the absolute magnitude of 
the most distant component My — 10. 4 m , which, ac- 
cording to the evolutionary tracks from Baraffe et al. 
(1997), corresponds to a mass of 0.22M Q for the metal- 
licity [Fe/H] = —2.0. Judging by the mass function from 
Chabrier (2003) derived for the Galactic halo stars, this 
mass is a mean for the halo field stars. Note that the dis- 
tance to the star used by Allen et al. (2000) is slightly 
larger, 180 pc instead of 170 pc that we adopted. We did 
not apply a correction to the mass estimate for this com- 
ponent due to the difference between the assumed helio- 
centric distances. 

Thus, denoting the primary component of the spec- 
troscopic pair, its closest companion, the speckle inter- 
ferometric component, and the most distant component 
with a common proper motion by A, B, C, and D, respec- 
tively (Fig. 2), we have M A w 0.67 M Q , M B w 0.24 M , 
M c w 0.33 M Q , M D w 0.22 M . 



5. ORBITAL PERIODS OF THE SUBSYSTEMS 
AND DYNAMICAL STABILITY 

Using the above mass estimates and Kepler's general- 
ized third law, we calculated the unknown periods of two 
subsystems of G89-14. The semimajor axis appearing in 
Kepler's law was determined for each subsystem from an 
empirical relation between the projected component sepa- 
ration and the orbital semimajor axis. Given the parallax 
of the system n and the projected separation between the 
components p, the expected semimajor axis can be calcu- 
lated from the formula (Allen et al. 2000) 



(a) = lo^+o- 146 , 



(2) 



where (a) is in astronomical units, p and tt are in arcsec- 
onds. 

Our calculations show that the speckle interferometric 
component and the SB1 pair make one revolution around 
the common center of mass approximately in 3000 yr 
(Pab-c), while the revolution period of the most distant 
component and the triple subsystem is about 650000 yr 
(Pabc-d)- For the outer subsystem, we took the expected 
semimajor axis (a) = 8565 AU from Allen et al. (2000). 
The ratio of the subsystem periods, 1 : 5769 : 1250000, 
is indicative of a high degree of hierarchy of G89-14 and, 
hence, its internal dynamical stability. Nevertheless, the 
outer subsystem has a low binding energy and is sub- 
jected to the destructive effect from giant molecular clouds 
and the stellar component of our Galaxy (Weinberg et al. 
1987). 
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Table 2. Initial conditions for constructing the Galactic orbit of G89-14 



V rad , km/s 


"/y r 


D, pc 


w, kpc 


2, kpc 


$, km/s 


W, km/s 


ft, km-kpc/s 


-36.9 


0.31 


170 


8.15 


0.03 


-165 


13 


194 



6. GALACTIC ORBIT 

To ascertain which subsystem of our Galaxy (halo, thick 
or thin disk) G89-14 belongs to, apart from the atmo- 
spheric metal abundance, it is also necessary to know the 
pattern of motion of the object in the Galactic gravita- 
tional field. The most complete description of the star's 
dynamical properties is the construction of its Galactic 
orbit. 

To construct the Galactic orbit, we must find a solution 
to the equations of motion in the gravitational potential 
of the Galaxy. For an axisymmetric potential, the spe- 
cific angular momentum is a conserved quantity and the 
motion can be described in a moving meridional plane. 
In cylindrical coordinates, the motion of the star in the 
meridional plane is described by a system of second-order 
differential equations: 



ft 
dt 



2 duj lu 3 ' 



(3) 
(4) 



d 2 z _ _d0 

dt 2 ~ 

At the same time, the motion of the meridional plane itself 
is specified by the equation 



dB_ 
~dt 



■—<) 1 



(5) 



where <p is the gravitational potential of our Galaxy, 
h = O • u> is the specific angular momentum of the body, 
w = VX 2 + Y 2 , 9 = arctg(T/X) and z = Z arc the cylin- 
drical Galactic coordinates, and O is the projection of the 
space velocity vector corresponding to coordinate 6; X, Y 
and Z arc the coordinates of the object in the rectangular 
coordinate system with the origin at the Galactic center 
and with the XY plane coincident with the Galactic plane. 
The second term on the right-hand side of Eq. (3) is the 
Coriolis acceleration acting on the body in the rotating 
frame of reference. 

The initial conditions in solving the equations of mo- 
tion (3)-(5) are the position and velocity of the star in the 
Galaxy. To specify the initial conditions, we must know six 
observable quantities: the object's coordinates (I and 6), 
proper motions in a (fj, a , " /yr ) and 8 (fig, "/yr), par- 
allax (w, "), and radial velocity (v rad , km/s). Table 2 
presents the input data for integrating the Galactic or- 
bit of G89-14. The projection of the space velocity vector 
onto to is denoted by <j>. In our calculations, we assumed 
the distance from the Galactic center to the Sun to be 
ujq — 8 kpc and the solar motion relative to the Local 
Standard of Rest to be (U®,V q ,Wq) = (-10.0, 5.0, 7.0 
km/s) (Dehnen and Binney 1998). Using formulas from 



Johnson and Soderblom (1987), we passed from the ob- 
servable quantities (ESA 1997; Latham et al. 2002) to the 
projections of the space velocities of G89-14 in the rect- 
angular coordinate system, (U,V,W) = (—165, —195, 13 
km/s). 

To calculate the parameters of the body's Galactic or- 
bit, it will suffice to solve two equations of motion, (3) and 
(4), numerically. For the solution, we used the seventh- 
order Runge-Kutta-Nystrom method with a variable step 
(see, e.g., Fehlberg 1972). The parameters of the derived 
Galactic orbit of G89-14 are presented in Table 3. As the 
gravitational potential of our Galaxy, we chose a three- 
component model described by Allen and Santillan (1991). 
Figure 3 shows the orbit of G89-14. Over 10 10 yr, the sys- 
tem has crossed the Galactic plane <~ 100 times. G89-14 
may have had a larger number of components that were 
lost through the tidal effect from the Galaxy, because part 
of the time the system moved in the Galactic plane, where 
this effect is maximal. We see from the figure that during 
its multiple revolutions, the star has closely approached 
the Galactic center (u> min = 0.3 kpc). The orbit is clearly 
chaotic and has a high eccentricity, which was determined 
from the formula 

Rapo Rperi 



R, 



apo 



R 



peri 



where R = \Juj 2 + z 2 ; R apo and R per i are the maximum 
and minimum distance of the star from the Galactic cen- 
ter. The orbital parameters, along with the low metallicity, 
allow us to confidently classify the object under study as 
belonging to the halo population. 

Carney et al. (1994) also obtained the parameters of 
the Galactic orbit for G89-14, suggesting that the star 
most likely belongs to the thin or thick disk. Given the low 
atmospheric metal abundance in the star, its belonging to 
the Galactic disk is less plausible. The main error in their 
calculations was the underestimated heliocentric distance 
(94 pc) determined by Carney et al. from the photometric 
parallax. 



7. CONCLUSIONS 

Our speckle observations at the 6-m BTA telescope re- 
vealed that G89-14 consists of four subdwarfs. A review 
of the literature showed that G89-14 with [m/H] = —1.9 
is the most metal-poor quadruple system known to date. 
This made it an interesting object for a more detailed 
study. 

Based on evolutionary tracks (Baraffe et al., 1997), 
available spectroscopic (Carney et al. 1994) and photo- 
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Table 3. Parameters of the Galactic orbit for G89-14 
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Fig. 3. Galactic orbit of G89-14. 



metric (Allen et al. 2000; Monet et al. 2003) data, and 
our speckle interferometry, we estimated the masses of the 
components of G89-14: M A w 0.67 M , M B w 0.24 M , 
M c ~ 0.33 Mq, M d « 0.22 M©. The ratio of the orbital 
periods of the subsystems Pab ■ Pab-c '■ Pabc-d = 
0.52 yr : 3000 yr : 650000 yr (1:5769:1,250,000) is indica- 
tive of a high degree of hierarchy of G89-14 and, conse- 
quently, its high internal dynamical stability. Using the 
three-component model of the gravitational potential by 
Allen and Santillan (1991), we constructed the Galactic 
orbit of the quadruple star under consideration (Fig. 3). 
Analysis of the Galactic orbital elements and the low 
metallicity of the system suggest that G89-14 belongs to 
the Galactic halo. Further observations of G89-14 by var- 
ious methods are needed to improve the orbital elements 
of the subsystems and the fundamental parameters of the 
system's components. 
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